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Summary. We describe a new model of massless thermal bosons which predicts 
an hyperbolic fluctuation spectrum at low frequencies. It is found that the partition 
function per mode is the Euler generating function for unrestricted partitions p(n). 
Thermodynamical quantities carry a strong arithmetical structure: they are given 
by series with Fourier coefficients equal to summatory functions at (n) of the power 
of divisors, with k = — 1 for the free energy, k = for the number of particles 
and k — 1 for the internal energy. Low frequency contributions are calculated us- 
ing Mellin transform methods. In particular the internal energy per mode diverges 
as = 2- with x = ^ in contrast to the Planck energy E = kT. The theory 
is applied to calculate corrections in black body radiation and in the Debye solid. 
Fractional energy fluctuations are found to show a 1/u power spectrum in the low 
frequency range. A satisfactory model of frequency fluctuations in a quartz crystal 
resonator follows. A sketch of the whole Ramanujan-Rademacher theory of parti- 
tions is reminded as well. 



1 Introduction 

According to the equipartition law of statistical mechanics, the available noise 
power P in the frequency interval dv is equal to kTdv pQ: this result is es- 
sentially Nyquist's theorem for the voltage noise (w 2 ) at a resistor R, i.e. 



P = -^r- = kTdv, where () means the average value. Since E = kT is 
the mean energy per mode, Nyquist proposed to add quantum corrections as 



energy. There are still controversies concerning the physical relevance of these 
relations: the Planck factor removed the ultaviolet divergence but this was 
reincorporated in the frame of quantum electrodynamics 

At the present stage, quantum statistical mechanics does not include in- 
frared corrections of the 1/u type. The infrared catastrophe was studied pre- 
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viously in non-stationary processes such as the scattering of electrons in an 
atomic field 0,0. 

Here we derive an alternative approach in which 1/v noise is a property of 
non degenerate bosons in equilibrium. We first observe that the quantum me- 
chanical partition function of a boson gas, with equally spaced energy levels, is 
the Euler generating function of p(n) : the number of indiscernible collections 
of the integer n. It relates to elliptic modular functions which are very exactly 
known. As a result the main contribution in the mean energy per mode is 
the infrared term ^ = instead of unity. The new law also leads to 1/v 
fractional energy fluctuations of the whole gas. 

Using the new approach and the density of states of the conventional 
approach we calculate the corrections to black-body radiation laws, including 
the density of photons, the emissivity and infrared fluctuations. We also apply 
the calculations to the phonon gas in a quartz resonator. 

The partition function Z of a non degenerate boson gas is given from 

lnZ=-]Tln[l-cxp(-/3e s )], (1) 

s 

where the summation is performed over all the states s of the assembly. In 
the conventional approach it is thus considered that the partition function Z 
per mode of frequency e s = hv s is such that \nZ = — ln[l — exp(— f3e s )]. 

In black-body radiation one accounts for the wave character of the quan- 
tum states by counting the number I of wavelenghts in a cubic box of size 
L 

l 2 = ^ = ll + ll + ll (2) 

where the summation ([Q) should be performed over all integers 1%, I2, h 
obeying @. 

This can be achieved by removing the discretness of energy levels and 
replacing the sum Q by an integral 

In Z = - / D{v) ln[l - exp(-/3/ii/)]di/, (3) 
Jo 

with D{u) = 2 x A7rV / the density of states [S]: the factor 2 happens due to 
the two degrees of freedom of polarization, c is the light velocity and V = L 3 
the volume of the cavity. 

From now we consider that to each mode is associated a set of equally 
spaced energy levels nhv, n integer, so that the partition per mode becomes 

InZ = - ^ ln[l - exp(-n/3H. ( 4 ) 

n>l 



As shown in Sect. (|2.1() this accounts for new multiparticle microstates not 
considered so far. 
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The summation above is well known in number theory and can be very 
accurately described using elliptic modular functions. At is will be shown, 
there are drastic consequences in the low frequency part of the spectrum, 
while the high frequency part is left unchanged. 

In the following the thermodynamical quantities will be defined as usual 



N = —-^TT, — t, the occupation number, (5) 

d(j3e a ) 

dlnZ 

= ~d{3~ ' internal energy, (6) 

dF 

S=—, the entropy, (7) 

din Z 

u = kT 2 — , the spectral energy density, (8) 

F=-kThiZ, the free energy, (9) 
dE 

e 2 = kT 2 ——, the fluctuations of the internal energy. (10) 



In all the paper the subscript ~ will indicate that we restrict the calculation 
to one single mode. 



2 Thermodynamics of the Euler Gas 
2.1 Euler generating function 

The partition function per mode Z in Q can be written in the Euler form |S] 

n>l J n>l 

with y — exp(— x) and x = pp. This is equivalent to the Boltzmann summa- 
tion 

z(y) = ^2 p ^ ex p(- nx )> ( 12 ) 

n>l 

where p(n) is the degeneracy parameter of the energy level nhv. It is known 
in number theory as the number of unrestricted partitions of the integer n, 
that is the number of different ways of calculating n as a sum of integers. 

For example with n = 4 we have p(4) = 5 and the corresponding indis- 
cernible collections are 

4 = 4 + + + (a), 

4=1 + 1 + 1 + 1 (b), 4 = 2 + 2 + + (c), 
4 = 3 + 1 + + (d), 4 = 2 + 1 + 1 + (e). 

(13) 
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This can be pictured in terms of the energy levels. The collection (a) means one 
particle on the level of index 4 and the three remaining particles on the ground 
state of index 0, i.e. Ahv = 1 x Ahv + 3 x Qhv. This collection is the only one 
considered in the conventional approach. The others microstates from (b) to 
(d) corresponds to different possibilities of bunching of the particles. Collection 

(b) means the four particles on the level of index 1, i.e. 4hv = 4x his, collection 

(c) means two particles on the level 2 and two particles on the ground state, 
i.e. Ahv — 2 x 2hv + 2 x Ohv, and so on. 

Properties of Euler generating function were studied in full details by Ra- 
manujan |7j in 1918 and completed by Rademacher [S] in 1973. An important 
result is the asymptotic formula 

p(n) ~ = exp(7r-y/2n/3) when n — > oo. (14) 

4nv3 

If instead of eq. (|llfl one uses the classical partition function 




Z — J2 n >i ex P(~ nhv/kT), one recovers the Planck's factor E/kT — x / (exp(x)— 
1). Such a dependence remains valid at high frequencies hvkT. If one intro- 
duces the zero point energy through the formula e Si „ — (n+ h)hi> with n > 1 
instead of e S! „ = nhv, then the low frequency part of the quantum spectrum is 
left unchanged but the ultraviolet catastrophe that one gets from the formula 
E/kT = x/(exp(x) — 1) + x/2 = a;coth(x/2), is cancelled as shown in Fig. 1. 
The new quantum theory thus solves simultaneously the questions concerning 
the 1/p power spectrum (which is observed but was not predicted) and the 
ultraviolet catastrophe (which is not observed but was predicted) jjj]. 
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2.2 Riemann zeta function and the free energy 
of the Euler gas 

The partition function Z(x) denned in Q is related to the Riemann zeta 
function £(s) through the Mellin transform as follows (Ref. [U3|, Eq. (6.3) p. 
464) 

/>OC 

r(s)C(s)((s + 1) = / (-InZix^x^dx. (15) 



(i 



Introducing <Jk(n) as the sum of k th powers of the divisors of n and the 
Dirichlct scries 

c(-)«*-*) = E f #. ( 16 ) 



n>l 



and computing the inverse Mellin transform one obtains (^H], p. 467) the free 
energy F as follows 

— — = — In Z(x) = ln(l — cxp(— nx)) = — <7_i (n) cxp(— ra). (17) 
n>l n>l 

with cr_i(n) = t7i(n)/n. There is thus a close relationship between the arith- 
metic of p(n) and that of divisors. This will be confirmed in the derivation of 
the others thcrmodynamical quantities. 

Since the main contributions to Z(x) are given by the poles at s = of 
r(s) and ((s + 1) and at s = 1 of C( s )> the free energy may be approximated 
in the low frequency part of the spectrum (^3> P- 58) 

F ir 2 1 , , x . x , „. 

with the error term X)z>i m (l — exp(— 47r 2 Z/x)). 

2.3 Dedekind eta function and the internal energy of the Euler gas 

One easily shows that the Mellin transform of the occupation number N(x) 
is r(s)((s) 2 . A similar derivation to the one performed in Sect. 12.21 leads to 

a m * — ' exp (nx) — 1 *■ — ' 

v ; n>l ^ y ' n>l 

In the low frequency region one gets (see p. 27) 

N(x) * Z^+J. (20) 
x 

where 7(2:) ~ 0.577 is Euler constant. 
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The Mellin transform of the internal energy E(x) is r(s)((s)((s — 1) and 
from the same method than above 



E[x) 



d{\nZ{x)) 



= hv E cM Z) -i = hv T. ^ eM-nx). (21) 

n>\ ^ v ' n>l 



An alternative derivation involves elliptic modular aspects. According to Nin- 
ham |1U) : All mathematics is a tautology, and all physics uses mathematics to 
look in different ways at a fundamental problem of philosophy - how to bridge 
the discrete and continuous. 

The link between the modular group and the Euler generating function is 
from the equality [E] 

^> = n = i>*" = 2E fP< w 

n>l y n>l ly ' 

where the domain of integration of Z(y) is taken to be the upper half complex 
plane of the new variable r 

y = exp(2Mrr), 9f(r) > 0. (23) 

Here we have 

~ / x x hv . , 

9f(r ) = — , x=—. 24 
2ir kl 

As shown in Sect. © Dedekind eta function acts on the full modular group 
SL(2, Z). At this stage we do not enter into the full ramifications of the theory 
and only emphasizes the connexion to the modular Eisenstein function 

G ^) = E? — ^< ( 25 ) 

' (tot + nY 

where the summation is performed over all non zero relative integers m and 
n and 3(r) > 0. It can be shown (^j, p. 29) that G^t) connects to the 
logarithmic derivative of rj(r) 

G 2 (r) = -4^ d(ln( ? (T))) , (26) 
dr 

with a Fourier expansion 

G 2 (t) = 2C(2) + 2{2mf ^ ax{n)exp{2mrvr). (27) 

n>l 

Using lO an d (EHJ- (E3> the relation ||3TJ is easily recovered. 
The low frequency expansion of internal energy is as follows 
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instead of the Planck result E ~ kT. 

One can also compute the entropy S from 

^ ^ j£ ^ 

- = — + \n(Z) = kT = o-i (n) (x + l/n) exp(-nx) 

n>l 

~ — + -\n— - - when x -> 0. (29) 

At very low frequency it results that the internal energy equals the opposite 
of free energy and one half the entropy. 



3 Application to Black-body Radiation 
3.1 Stefan-Boltzmann constant revisited 

The Stefan-Boltzmann constant is an integrated measure of the emissivity 
of a black body In the conventional approach the partition function is 
calculated from the integral using the density of states D(v) = 87T ^ U that 
is 

\nZ = 8%V\^j x 2C(4), (30) 
with C(4) = 7r 4 /90 and we used the Mellin integral formula 

C(s + 1) = - 7T -/ x s - 1 ln(l-exp(-x))dx. (31) 

r{s) Jo 

The Stefan-Boltzmann constant erge is defined from the free energy 

F = -kT\nZ=-^VT i witha SB = ^^. (32) 
3c 15c z h J 

If instead of (|30|l one uses the general formula 



ST*** 



InZ = -8ttV — / x z > ln(l - exp{-nx))dx, (33) 



the interchange of the integral and the sum leads to 

\nZ = -8nv(^-\ ,2,(4) I > ' — I ■ CUl 



ch 



As a result we find a free energy (and a modified Stefan-Boltzmann constant) 
in excess with a factor $Z n>1 -tf = C(3) — 1.20. If one uses the alternative 
derivation in terms of the divisors one recovers the mathematical formula 1161) 
with s = 3 and k = — 1. 
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3.2 The density of photons 

The number of photons in the bandwidth dv is 

dN(v, T) = D(v)N{ Vl T)dv, (35) 

with the occupation number N(v, T) — (exp(f3hu) — 1) _1 in the conventional 
approach. Integrating one gets per unit volume 

n( t Ks _ 8 j^)\ 2m (36) 



V \ ch J jq exp(x) — 1 \ch 

In the general approach the occupation number is defined from the sum- 
mation (ftT?|) and we need to evaluate 



x 2 



E^^=2C(3)(|En^=2C(3f. (37) 



This is in excess of a factor £(3) ~ 1.20 as for the free energy. If one compares 
the calculation in terms of divisors one recovers the mathematical formula 
JHJ with s = 3 and k = 0. 

3.3 Planck's radiation formula revisited 

The energy within the bandwidth dv is defined as 

dE(v, T) = D{v)E{ Vl T)dv = u(v, T)dv, (38) 

with E{y,T) = hv{eyL~p{l3hv) — 1) _1 in the conventional approach and with 
u(v, T) the energy spectral density. We get the Planck's radiation formula 



8nhV 

u{v, T) = — — — ■ — -. 39) 

c J exp{phv) - 1 



The black-body emissivity is defined as 

e b {v,T) = -^u{v,T). (40) 

At very low frequency the conventional result leads to the Rayleigh- Jeans 
formula 

[e b {v,T)] KJ = 2^v 2 T, (41) 

which is independent of Planck's constant and is proportional to the inverse 
of the square of wavelength A = c/v. 
In the new approach the emissivity is 
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e b (u,T) = — V 42 

>i ex P\ n P' lv ) ~ 1 

At very low frequency one uses (|28[) with the result 

M,,T)] LF = ^4^T 2 . (43) 

Thus the v 2 T dependance is replaced by the vT 2 dependance, the low fre- 
quency emissivity now depends on the Planck constant and on the inverse 
wavelength; there is a ratio between the new result and the one predicted 
by the Rayleigh- Jeans formula. 



3.4 Radiative atomic transitions 



Let us now consider the equilibrium between atoms and a radiation field, 
allowing the emission or absorption of photons of frequency 

v = V\ — U2, (44) 

where ti — hv% is the energy in the upper state and e% = hv\ in the lower 
state. 

The conventional theory, as derived for the first time by Einstein, states 
that the rate at which atoms make a transition 1 — > 2 in which one photon is 
absorbed is equal to the rate at which atoms emit photons, so that 

B l2 N lU (y) = A 2 iN 2 + B 2 iN 2 u(v), (45) 

where N\ and N2 are the occupation numbers of atoms in levels 1 and 2, A21 
is the spontaneous absorption rate, B12 is the induced emission rate and u(v) 
is the energy density in the radiation field as given in (|39|l . 

In thermal equilibrium the occupation numbers in states 1 and 2 obey the 
Boltzmann law 

^=exp(--). (46) 
Using (|46|l and i|39|) one gets the well known formulas 

A21 = A, B X 2 = B 2 i=B and ^ = , (47) 

where A = cjv is the wavelength of the radiation field. 

If one uses the general formula one gets low frequency corrections in the 
spontaneous to stimulated emission ratio. Using the low frequency expression 
l|28|) for the internal energy this yields 
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The A/B low frequency ratio now depends on the inverse of the square of 
wavelength A and is independant on the Planck constant h, in contrast to 
the h/X 3 dependance of the standard ratio. The spontaneous to stimulated 

2 

absorption rate is enhanced a factor 5- over the conventional one. 



3.5 Einstein's fluctuation law revisited 



According to the conventional Einstein's approach |13j the energy fluctuations 
of a system in equilibrium within a larger system of temperature T are 



((E-(E)) 2 ) = - 



d(E) 
d(3 



k,T 



dT ' 



(49) 



One can reformulate this relation for the fluctuations of the energy dE 
u{v,T)dv in the bandwidth dv 



de 2 = kT 2 —dv = S u {v)dv, 



(50) 



with u the energy spectral density and S u (v) the power spectral density of 
the fluctuations of u. One gets (^|, p. 429 ) 



de 2 = {hvu - 



iirv 2 V 



u 2 )Vdv. 



(51) 



The first term at the right hand side is the one corresponding to the high fre- 
quency part of the spectrum (Wien's law): it is of a pure quantum nature and 
is corpuscular like; the second one corresponds to the low frequency (Raylcigh- 
Jeans) region: it is purely classical and wavelike. In the low frequency part of 
the spectrum they are fractional energy fluctuations of the random walk type 



1 



R.7 



8ttV v 1 



(52) 



In the new approach one uses the low frequency energy density u(v, T) 



so that instead of (g>2J) one gets 



S u {v) 



LF 



3 he 3 1 
2^n?fc7V' 



(53) 



This is the announced quantum \ jv fluctuation spectrum. There is a reduced 
low frequency noise and the ratio between the new result and the Einstein- 
Rayleigh- Jeans one is i^r- 
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4 Application to a Phonon Gas and to the 1/ f 
Frequency Noise of a Quartz Resonator 

4.1 The specific heat of a phonon gas revisited 

The properties of the phonon gas are quite similar to those of the photon 
gas except for the new form of the density of modes as g(v) = 12 $ v v 2 with 

-§- = -g- + ^ where c p h represents the average wave velocity and c t and 

ci are the transverse and longitudinal velocities for an isotropic solid [Jj]. The 
maximal vibrational frequency v m (Debye frequency) is defined from the total 
number of allowed quantum states 



3N = / i/dv that is v m = — — , (54) 

C ph Jo \ 47rV J 

where Nq is the number of atoms in the volume V. 
In the conventional theory [d] we get 

In Z ph = ^ / v 2 In [1 - exp(-/3H] dv. (55) 

v m Jo 

The internal energy follows from the formula 

9RT 

E=^D(x m ), (56) 
x m 

and the constant volume specific heat C v = dE/dT equals 

C v = 3R [D(x m ) - x m D'(x m )} , (57) 

with D(x m ) = Jo™ c ^)-i the Deb y e function, and x m = ^, with 6> D = 
hu m /k the Debye characteristic temperature. 

The case T9d, D(9t>) — > 1 corresponds to the Dulong-Pctit value C v ~ 3R. 
At very low temperatures one gets the cubic temperature dependence C v ~ 

¥x M (x) 3 . 

In the new approach 

In Z ph = ^ / Vv 2 In [1 - exp(-n/3hu)] dp. (58) 

1/0 «>i 

Debye results are found unchanged except for an extra multiplicative factor 
in the specific heat as was the case for the integrated emissivity in Sect. 13.11 
that is 

Si™ _ £( 3 ) „ i.20. (59) 

At very low temperatures the electronic contribution to the specific heat which 
decreases as T, dominates the lattice contribution, which decreases as T 3 . This 
is accounted for in the conventional way. 
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4.2 1/f noise in a quartz resonator 

Specific heat is involved in the energy fluctuations of a canonical ensemble 
from the relation 

e 2 = kT 2 C v . (60) 

2 / , \i/2 

The relative energy fluctuations follows as -gy = ( j^- \ which is on the 

order 10 -11 for iVo = 10 23 , the Avogadro number. 

For energy fluctuations in the bandwidth dv, the main difference with the 
conventional theory lies in the low frequency region, as was the case of the 
photon gas. We find the quantum 1/v formula 



S u (y) 



= ^J_ = A* (61) 
LF 4tt 3 V kTv Vv K > 

The method can be used to predict fractional frequency fluctuations in a 
quartz crystal resonator from the formula 1 |14 | -|15 | 

w 2 4Q 4 Vv v ' [ ' 

where to and Q are the frequency and quality factor of the resonator. Using 
c ph ~ 3.5 x 10 3 m/s, we find A ph ~ 5 x 1CT 4 . For a 5 MHz P5 quartz 
crystal resonator with Q ~ 2 x 10 6 , the active region under the electrodes has 
thickness t = 5A/2 ~ 3 mm, and section S ~ 3 cm 2 , that is V ~ 1 cm 3 . The 
resulting 1/v factor is h-i = A q4 V ~ ~ 6 x 10 -24 . This is the order of 
magnitude found in experiments |14j . 



5 Ramanujan— Rademacher Theory of Partitions: 
a Short Reminder 

Besides the low frequency approximations encoutered in Sect. @ there is a 
an exact method to calculate the number of partitions p(n) first discovered 
by Ramanujan and improved by Rademacher |Hj thanks to an integration 
along Ford circles in the complex half plane. For completness we remind here 
the main points of the theory from which the results in Sect. (J2J may also be 
derived . 

From well-known mathematical arguments [7] (p. 113) one can get the 
leading term for the case < y < 1 and y — ► 1 from the expression 2 . 



1 To establish the formula one writes the equation for a lossy harmonic oscillator 
and one postulates that the 1/v fluctuations are present in the loss coefficient. 

2 We have 
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The use of y = exp(— hv/kT) corresponds to the low frequency approximation 
at v — ► 0, that is l — y= 1 — cxp(— hv/kT) ~ hv/kT. This leads to the leading 
low frequency term in the free energy (|18|1 and internal energy H2KJI. 

They are similar formulas associated with rational points which are located 

at 

P 

y pq = exp(2?7r-), (64) 

on the unit circle \y\ = 1. The leading term in the expansion of Z(y) corre- 
sponds to the fundamental mode | = j. 

The general method to compute rational contributions is a master piece 
of twentieth century mathematics (|Zj), (0)- It uses the connexion of Z(y) to 
the elliptic modular functions. 



5.1 The fundamental contribution 



To compute the contribution of the fundamental point 1/1 of the unit circle 
\y\ = 1 one uses the property (H§|, p. 96, (TT], P. 58) 



z(y) 



V ,(lni) 1 / 2 exp 



'2it 



y 



61n- 



Z(y') with y = exp 



4n 2 



lny 



(65) 



In the low frequency region y — exp(— x) ~ 1 so that y' = exp(— 4ir 2 /x) ~ 
and Z(y') ~ 1. Low frequency approximations of the free energy (|18|l and of 
the internal energy (|28|l follow. There are similar formulas associated with the 
other rational points of the circle as shown below. 

To get the leading term in p(n) one uses the Cauchy formula 



If Z(y) , 
P(n) = 77- f -T^dy, 

where <f means an arbitrary closed loop encircling the origin 

iogz (y) =5>g-^=£^=]r— 

i— i 1 — y n t-^i m ^-^ mil 



(l-ymy 



and 



so that 



my 



\l-y)<l~y m <m(l-y), 



1—1/ m^ 1 — y L — ' 



y 



y <-^t m- L-y 

m in 

Each of the above series has the limit 7r 2 /6 when y — * 1 and so log Z(y) 



(66) 



6(1-1/) 
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Substituting l|65|) in l|66[l with Z(y') = 1 one can obtain 



This includes l|14|l in the limit n — > oo but is much more accurate. 
5.2 Farey contributions and Ford circles 

From now on we extend the domain of definition of Z(y) to the complex plane 
and we introduce the new variable r and Dedekind eta function rj(r) as defined 
in 

It can be shown that j](t) is a modular form of degree —1/2 on the full 
modular group. It acts on the generators of such a group through the relations 

EH 3 

r)(r + 1) = cxp(i7r/12)77(r); f/(-l/r) = (»/A) 1/2 ?7(t). (68) 
To express the partition function one uses the Cauchy formula 4 

1 /" Z(v) f~o+l 

pin) = <f> — Wdy = / Z [exp(2j7rr)] exp(— 2iTTTn)dT. (69) 

2^7^ / J ro 

In the third term above this corresponds to a path of unit length starting at 
an arbitrary point in Ti. 

The choice of the integration path comes along in a natural way by using 
the connexion of Z(y) to the modular group. Let us observe that the set of 
images of the line r = X + i, X real, under all modular transformations 

T = -. with p,p',q, q' integers and |pq' — q'p| = 1, (70) 

qr + q 

can be written as 



3 For more general modular transformations, we have 



,pr+p , , qr + q' , , 

^qVTq 7 ,Q \l i 

with e a 24t/i root of unity related to Dedekind sums as defined in 17911 . See 
p. 160. 

4 Let f(z) be an holomorphic function of the complex variable z. The Cauchy 
formula for the derivatives is as follows 

r(a) = — / , f{Z )+, dz 

where <f> means a closed contour encircling the pole. It is applied in Ho 91 to the 
partition function / = Z with a — so that J—Lpl — p(n) 
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( P - + — 



1 

V' 



2«- 



(71) 



with radius 



Equation defines circles C(p, q) centered at points r = ^ 
l/2q 2 . They are named after L.R. Ford who first studied their properties in 
1938 jS]. Ford circles are easily generated by using the ordered Farey sequence 







< 



< P_l < P_L 



('li 



<h + 92 



_P2 P2 



< 



'12 



(72) 



To each - belongs a Ford circle in the upper half plane, which is tangent 
to the real axis at t = |. It can be observed that Ford circles never intersect. 
They are tangent to each other if and only if they belong to fractions which 
are adjacent in some Farey sequence. 



If 



< E < £2. are three adjacent fractions in a Farey sequence then 



C(p,q) touches C(pi,qi) and C(p2,<?2) respectively at the points 



t L -P+fL and t R - H + 

. „„ - + (, pq and r pq - + C, pq , 



pi 



(73) 



where 



5i 



q(q 2 + ql) q 2 + q\ 



and C^q = 



C12 



q(q 2 +q|) 



q 2 + ql ' 



(74) 




i+1 



q 1/4 1/3 1/2 2/3 3/4 ^ 

Fig. 2. Rademacher's path of integration. 



In Rademacher's approach (which improves Ramanujan's one) the unit 
length path on Ti. is choosen so as to go along Ford circles 
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p(n) 



(p,g)=i 

E 

0<p<q<N J 1-Pl 



Z [exp(2i7rr)] exp(— 2iTTTn)dT, 



(75) 



where 7 pq is the upper arc on a Ford circle which connects points of tangency 
at r^ q and t*. 

Each Ford circle C(p,q) is parametrized by the expression r = | + £, where 

the variable £ runs on an arc of the circle l^-^rl = ^2. If one uses z such 

that C = *£) a Ford circle is mapped onto the circle \z — || = | and lf75|) 
transforms as 



p(ra) 



(p.g)=i 

{— exp(-2z7rn-) 



0<p<g<JV 



exp(27ri 



ex P( — 



where and z^ Q follows from l(75)l. 



(76) 



5.3 Farey contributions and Dedekind sums 

To compute (|76|) one uses the transformation formula 11' 21) . After some manip- 
ulations and using z/q instead of z(see Ref. jS], p. 269), one gets the formula 
which generalizes (|65f) 



1/2 

Z(y)=uj pq [-\ exp(-- 



12z 12q 2 



)Z(y') 



with 



2inp 2ttz . 2mp' 27r . 

y = exp( — ), y = exp( )and pp 

q q z q z 



(77) 



-l(mod q). (78) 



The so-called Dedekind sums s(p, q) are introduced by 

q 



= exp (iTTs(p, q)) with s(p, q) = ^ 



1=1 



pi 

q 



(79) 



where [ ] in l|79|l denotes the integer part. 

For the calculation of p(n) one uses an approximation similar to the one 
used in l|65|) . If z is a small positive real number, then y is near exp(2i7r£), 

the modulus at that point \y'\ = exp(— — ) ~ and Z{y') ~ 1. 
As a result (|76|) can be readily integrated and the final result is 

1 



q>l 



Ar, 



(80) 
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withK q = ^| A n = ^/n-i_, 

and A q (n) = Lu pq exp(— 2i7rn— ). (81) 

p mod(q) ~ 



6 Conclusion 

We have found that quantum 1// noise may be a property of photons or 
phonons in thermal equilibrium. The theory connects elliptic modular func- 
tions and quantum statistical mechanics as in superstring theory, but with 
physical relevance in the macroscopic realm of infrared divergences of solid 
state physics. Main results are an enhanced energy per mode at low frequency 
(one gets ^ ~ \hJ instead of the Planck result E ~ kT) and an enhanced 
integrated radiation, photon density and phonon specific heat (with a factor 
of C(3) — 1.20). Low frequency fluctuations are reduced and one obtains frac- 
tional energy fluctuations with al/i/ spectrum in contrast to the random walk 
1/v 2 of the standard theory. 

Main steps of the general mathematical theory of p(n) is reminded in 
the last section. It is based on the fact that there are rational singularities 
on the unit circle in the partition function. The fundamental mode in the 
Farey decomposition leads to a satisfactory account of the infrared part of the 
spectrum and remains correct in the high frequency region up to x = x\ = 4tt 2 . 
The physical meaning of higher modes has still to be understood. 

Farey series and Ford circles were used recently in the different context 
of l/f noise in phase locked loops |TJ|. In this last case they connect to the 
arithmetical functions found in prime number theory. 
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